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ABSTRACT 

We present a survey of ~ 800 square degrees of the galactic plane observed with the QUaD telescope. The 
primary product of the survey are maps of Stokes /, Q and U parameters at 100 and 150 GHz, with spatial 
resolution 5 and 3.5 arcminutes respectively. Two regions are covered, spanning approximately 245-295° 
and 315-5° in galactic longitude /, and -4 < b < +4° in galactic latitude b. At 0.02° square pixel size, the 
median sensitivity is 74 and 107 kJy/sr at 100 GHz and 150 GHz respectively in /, and 98 and 120 kJy/sr for 
Q and U. In total intensity, we find an average spectral index of a = 2.35 ± O.Ol(stat) ± 0.02(sys) for \b\ < 1°, 
indicative of emission components other than thermal dust. A comparison to published dust, synchrotron 
and free-free models implies an excess of emission in the 100 GHz QUaD band, while better agreement is 
found at 150 GHz. A smaller excess is observed when comparing QUaD 100 GHz data to WMAP 5-year 
W band; in this case the excess is likely due to the wider bandwidth of QUaD. Combining the QUaD and 
WMAP data, a two-component spectral fit to the inner galactic plane {\b\ < 1°) yields mean spectral indices 
of as = -0.32 ± 0.03 and = 2.84 ± 0.03; the former is interpreted as a combination of the spectral indices 
of synchrotron, free-free and dust, while the second is attributed largely to the thermal dust continuum. In 
the same galactic latitude range, the polarization data show a high degree of alignment perpendicular to the 
expected galactic magnetic field direction, and exhibit mean polarization fraction 1 .38 ± 0.08(stat) ±0.1 (sys)% 
at 100 GHz and 1.70±0.06(stat)±0.1(sys)% at 150 GHz. We find agreement in polarization fraction between 
QUaD 100 GHz and WMAP W band, the latter giving 1 . 1 ± 0.4%. 

Subject headings: Surveys — Submillimeter — cosmic microwave background polarization — cosmology: 
observations — diffuse radiation — Galaxy: structure — Galaxy: disk — ISM: sti^ucture 



1. INTRODUCTION 

Radio and sub-mm observations of the galactic plane yield 
important insights into many astrophysical processes associ- 
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ated with galaxies like our own, from the large scale proper- 
ties of magnetic fields, to smaller scale phenomena associated 
with star formation. 

These properties can often be inferred from low resolution 
observations of diffuse galactic components, or large samples 
of representative objects distributed through the galaxy. Three 
mechanisms contribute to the diffuse galactic emission in to- 
tal intensity in the radio and sub-mm: Synchrotron radiation, 
which dominates below ^ 60 GHz and is caused by relativis- 
tic electrons spiralling in magnetic fields; free-free emission, 
which is generated by non-relativistic electron-ion interac- 
tions; and radiation from vibrational modes of thermal dust, 
whose emission dominates above ^ 100 GHz. Of these, syn- 
chrotron and dust are appreciably polarized, with typical po- 
larization fractions close to the galactic plane of ^ 2-4% and 
~ 1 % respectively ( |Kogut et al.|2007| l; both result in polarized 
light aligned perpendicular to the magnetic field. 

The polarization of dust is due to prolate grains align- 
ing with their long axis perpendicular to the local magnetic 
field (L azarian||2003[ ), with the polarization fraction depen- 
dent on the grain size distribution and their overall align- 
ment (e.g. Prunet et al. 1998 ). Observations of polarized 
starlight via dust absorption have indirectly demonstrated a 
large degree of coherence of the magnetic field in our galaxy 
and others ( Heiles 1996 ; Zweibel & Heiles|1997[ ). However, 
these measurements can be biased by lines of sight with low 
column densities ( |Benoit et al.|2004l ) and therefore more di- 
rect probes of the galactic dust are desirable, not only for the 
study of dust in its own right, but as a probe of the magnetic 
field responsibl e for most large-scale p olarized emission in 
the galaxy (e.g. Hildebrand et al.|[2000| l. Sub-mm polariza- 
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tion vectors can be reasonable tracers of the magnetic field 
structure even for relatively dense clouds, and are therefore 
an attractive option for this line of study. 

The emissive properties of dust in the sub-mm have also 
attracted attention from the CMB community, since diffuse 
galactic polarization poses a challenging obstacle to the detec- 
tion of primordial gravitational waves via the 'B-mode' polar- 
ization signal (e.g. ' Hu & White][T997l [Dunkley et aT][2009l l. 
As the B-mode power spectrum is predicted to peak on an- 
gular scales ^1°, characterization of dust as a foreground is 
essential to account for this 'contaminant' from surveys over 
larg e areas of sky, suc h as that expected from the Planck satel- 
lite (|WiaetaI]2002|. 

Observational constraints on diffuse galactic polarization 
from dust are currently limited to a small number of experi- 
ments, including WMAP from 23-94 GHz at resolution up to 
- 0.25° ( |Kogut et al .120071 [Gold et al.|2009|l, and Archeops 
at 353 GHz smoo thed to ~ 1" resolution ( |Benoit et al.|2004[ 
|Ponthieu et al.|20 05 ). While a variety of models exist for es- 
timating the unpolarized contribution of dust (e.g. Finkbeiner' 
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et al. 1999, hereafter FDS), the limited number of experiments 



at dust-dominated frequencies has prevented detailed compar- 
ison to observations. Furthermore, the lack of angular reso- 
lution of such experiments means emission from diffuse and 
discrete sources cannot be separated, particularly in the plane 
of the galaxy (|/7| < 10°). 

The study of discrete sources at and above 100 GHz yields 
insights into the process of star formation. In star forming 
regions, thermal dust efficiently absorbs UV light from star 
formation, and re-radiates in the sub-mm where dust is op- 
tically thin. Observations near the spectral peak can there- 
fore probe the centers of dense cores and constrain the stellar 
core mass function (e.g . [Netterfield et al.|2009t|Schuller et al.| 
[2009i|Olmi etal.|2009D . 

In this paper we report an ^ 800 square degree survey of 
the galactic plane with the QUaD telescope, which operated 
at 100 and 150 GHz with angular resolution of 5' and 3.5' re- 
spectively, in Stokes /, Q and U parameters. A survey of this 
size, frequency and angular resolution can be used to investi- 
gate the polarized and unpolarized properties of both diffuse 
emission and discrete sources. The mapmaking and properties 
of the diffuse emission form the core of this paper; a compan- 
ion publication (hereafter the "Source Paper" — Culverhouse 
et al. in prep.) contains analysis of the compact source distri- 
bution in the survey. 

2. INSTRUMENT SUMMARY AND OBSERVATIONS 

Here we summarize the features of the QUaD experi- 



ment — a detailed description can be found in (Hinderks 
et al.|2009] l, hereafter referred to as the "Instrument Paper". 
QUaD was a 2.6 m Cassegrain radio telescope on the mount 
originally constructed for the DASI experiment (Leitch et al. 
2002] !, and enclosed in an extended reflective ground shield. 
The receiver consisted of 31 pairs of polarization sensitive 
bolometers or PSBs ( Jones et al.|2003| ), 12 at 100 GHz, and 19 
at 150 GHz, with each PSB pair located within a single feed. 
The PSB pairs were split into two orientation groups to allow 
simultaneous measurement of Stokes Q and U. QUaD oper- 
ated from February 2005 to November 2007; the observations 
reported in this paper were taken over 40 days between July 

and October 2007 . 

As with the QUaD CMB observations ( Ade et al.|[2008l 




Pryke et al.|, 2009^ Brown et al.^ 2009) , the second of which 



we refer to as "P09", a lead-trail field differencing scheme 



Fig. I. — Top.' Coadded PSB pair difference map of a single 150 GHz feed 
from one day of observations, with no field difference. The nine rectangular 
lead-trail pairs of fields observed on this day are shown in blue and red re- 
spectively. Note the strongly polarized ground pickup which repeats between 
lead-trail pairs. The colour scale is MJy/sr, with the green line indicating the 
plane of the galaxy. Green squares indicate the starting point of the scan pat- 
tern for each field. Bottom: As above, but with field differencing applied — 
the ground signal is now heavily suppressed. 

was employed to allow subtraction of contaminating ground 
pickup. Each day, nine lead-trail pairs of fields were observed, 
with the lead field tracking center aligned with the plane of the 
galaxy b = 0. For each lead field, the companion trail field re- 
peated exactly the same azimuth (az) and elevation (el) scan 
pattern with respect to the ground, but with the tracking cen- 
ter Ihr later in RA. In order to minimize the possibility of 
temporal variation in ground signal, each trail field was ob- 
served immediately after its companion lead field, resulting 
in a lead 1, trail 1, lead 2, trail 2 etc ordering — see Figure [T] 
for a plot of this scheme. Though the signal levels in trail 
fields are generally much smaller than the lead field, some 
spurious signal is introduced into the lead field due to field- 
differencing. A discussion of this effect is presented in Ap- 
pendix |A3] while Appendices |B . 1 1 and |B .2| demonstrate that 
the recovery of global parameters of the diffuse emission are 
unaffected by field-differencing or filtering. 

Whilst tracking each field center, the telescope scanned for- 
ward and backward in az with a total az range of 15° at a rate 
of 0.4°/s, followed by a step of -0.02° in elevation. This was 
repeated 35 times to build up a raster map of the sky, before 
slewing to the next field where the same scanning pattern was 
repeated. Each field took approximately Ihr observing time, 
and covered 0.7° in dec, a total of 6.3° per day. Figure [T] 
shows a graphic representation of the scanning strategy, and 
demonstrates that ground signal is cleanly removed by field- 
differencing. 

The entire survey coverage is shown in Figure [2] and was 
limited in declination range by two factors. First, the beams 
intersect the ground shield at elevations lower than ^ 25° 
(dec > -25°). Second, at dec < -60° the galactic plane is 
nearly parallel to the horizon — filtering the scans, neces- 
sary to remove atmospheric contamination, would remove the 
bulk of the diffuse galactic emission unless the scan length 
was massively increased, resulting in a loss of sensitivity. 

Two portions of the galaxy are available between -60 
and -25° in dec, between 110° < RA < 190° and 210° < 
RA < 290°. In galactic longitude and latitude {l,b) these cor- 
respond to a latitude range approximately -4° < b < 4°, and 
240° < / < 295° and 315° < / < 5° — hereafter, these two 
regions are loosely termed the 'third quadrant' and 'fourth 
quadrant' . With our coordinate constraints, the third and fouth 
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Fig. 2. — Total QUaD galactic survey area at 150 GHz (black lines) plotted over WMAP Syr W band total intensity map )Hinshaw et al.|2009) . The color scale 
is in MJy/sr. For one day of data, lead field coordinates for center pixel are plotted in dark gray; trail fields are in light gray. Ihis lead-trail observing strategy is 
repeated over the rest of the survey region within the black locus. 



quadrant regions each spanned 32° in dec; the 40 days of ob- 
servations were equally divided between the two regions, al- 
lowing four complete passes over each. 

3. LOW LEVEL DATA PROCESSING 

Low-level processing of the raw timestream is performed 
using the same steps described in P09. The timestream is 
deconvolved to remove the effect of the bolometer timecon- 
stants, deglitched, and downsampled. Between observations 
of each field, relative gains between FSB pairs in a feed, and 
between feeds in each frequency group, are determined us- 
ing 'el-nods'. In this calibration procedure, the telescope is 
moved first up then down by one degree in elevation — this 
injects a large signal into the timestream due to the common- 
mode atmospheric gradient; further details are given in the 
Instrument Paper. 

As seen in Figure [T| the ground pickup is strong compared 
to the sky signal of interest. A misalignment between abso- 
lute lead and trail field azimuth coordinates Saz causes ground 
signal to cancel imperfectly, and can lead to significant con- 
tamination for (5az > 0. 1 ° . This is carefully corrected by using 
the pointing information to realign the lead and trail fields on 
a scan-by-scan basis, at the expense of a small quantity of data 
where there is no overlap between lead and trail scans. Ap- 
proximately one out of nine lead-trail pairs require correction 
of up to 0.4°, leading to a loss of data of order 0.1%. 

Further data is rejected from visual inspection of field- 
differenced maps, made using sum and difference data from 
each feed. Maps of the same fields taken on different days are 
compared to distinguish the repeatable sky signal from spu- 
rious contamination. The pair difference maps in particular 
are useful because the amplitude of the contamination, which 
appears as spurious 1// noise, is considerably larger than that 
of the polarized galactic signal. Typical rejection rates are 
one out of nine fields for four bolometer pairs, a loss of data 
of ~ 0.75%. 

4. MAPMAKING 

The map-making scheme used here is a multi-stage adapta- 
tion of the 'naive' mapping used in P09, and requires infor- 
mation on telescope pointing (both absolute and the relative 
offsets of each PSB in the focal plane), and the PSB angles 



and efficiencies to construct the Stokes /, Q and U maps. As in 
P09 and unless stated otherwise, the pixelization for all maps 
in this paper is in RA and dec, using 0.02° square pixels. Fur- 
ther details of constructing /, Q and U maps from timestream 
may be found in P09; here we summarize the basic points. 

4.1. Pointing, PSB Angles and Efficiencies 

Absolute pointing was determined from a nine parameter 
online pointing model, derived from optical and radio ob- 
servations as described in the Instrument Paper. This was 
shown to have an absolute accuracy over the hemisphere of 
^ 0.5' rms from pointing checks on RCW38 and other galac- 
tic sources, taken over two seasons of CMB observations. 

The scatter in the centroid positions for a given detector rel- 
ative to the boresight was consistent with the overall pointing 
wander of ^ 0.5' rms, and the offset angles of each detec- 
tor around the focal plane showed no evidence for systematic 
changes with time. The detector offset angles used in map- 
making are the mean of the values observed from RCW38, 
with an estimated uncertainty of ^ 0.15'. 

PSB polarization angles and efficiencies were determined 
using a chopped thermal source placed behind a polarizing 
grid and observed at many angles — further details are given 
in P09 and the Instrument Paper. From these observations we 
also measure our mean cross polar leakage (the response of a 
single PSB to anti-aligned radiation) as e = 0.08 ±0.015. This 
mean value is applied to all channels when constructing maps, 
and all simulations include the scatter about the mean. 

For an experiment of this type cross polar leakage does not 
imply leakage from total to polarized intensity — it is simply 
a small loss of efficiency, which is corrected by an additional 
calibration factor applied to the polarization data. 

4.2. Mapmaking Algorithm 

Before binning into maps, the data is first field-differenced 
to remove ground contamination — this operation is per- 
formed directly in the timestream. For each feed the sum and 
difference of the data is then taken for each pair of PSBs. To 
construct the / maps, we coadd the pair-sum data for each 
feed on each day. For polarization, a 2 x 2 matrix inversion 
is required for each map pixel to convert from the pair dif- 
ference data to Q and U maps. This matrix expresses the po- 
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larized sky intensity projected onto each PSB, which is mea- 
sured in the pair-difference timestream. To invert such a ma- 
trix we require each pixel be measured at two PSB angles — 
this is achieved with the two orientation groups in the QUaD 
focal plane. However, before coadding the data into maps, 
timestream filtering is required to reduce the low-frequency 
noise, which causes striping in the scan direction in the maps. 

4.2. 1 . Initial Filtering 

In addition to detector noise, bolometer drifts and atmo- 
spheric l/f noise are a large contribution to the sum data 
(since the atmosphere is largely unpolarized, its contribution 
is common-mode and is thus heavily suppressed in the differ- 
ence data). The CMB analysis of P09 subtract a third-order 
polynomial from the timestream to limit the effects of this 
l/f noise; though this filtering removes sky signal, the effect 
is accounted for in simulations, which is a feasible method for 
a power spectrum analysis. However, we are interested in the 
spatial distribution of the galactic signal, and the timestream 
may not simply be filtered in the same way. For example, in 
the sum data the bright galactic emission will dominate the 
polynomial fit, leading to regions of unphysical negative sig- 
nal when the polynomial is subtracted from the data; a min- 
imal level of filtering is therefore desirable. Here, the end 
portions of each scan — the most distant parts of the scan 
from the galactic plane — are used to determine a DC level 
and slope, which is then subtracted over the entire scan. The 
same procedure is used for total and polarized intensity data. 
This choice of filtering scheme effectively forces the maps to 
be zero at the edges, a consequence of our inability to measure 
the DC level of the sky brightness. 

The amount of scan ends used is a trade off between better 
determination of the 1 // noise, which requires an increasing 
fraction of the scan, and larger regions of negative intensity in 
the final maps. The second of these effects arises because the 
DC level subtracted from the timestream is influenced more 
by the bright galactic signal as more of the scan is used. 

After some experimentation, using a scan fraction /, = 25% 
(i.e. ^ 12.5% at either end) was found to result in total in- 
tensity maps with few negative regions, while keeping resid- 
ual 1 // noise to a minimum. The quantitative results on the 
diffuse emission, presented in Section|5] are not significantly 
changed by adopting a different fraction of the scan, or using 
a mask of fixed width either side of the galactic plane. 

From the sum/difference data we construct a minimally- 
filtered map, termed the 'initial map' or mo: scans are filtered 
by removing a DC-level and linear slope, as described above. 
The timestream is then coadded into the map using the point- 
ing information for each pair of PSBs, and weighted by the 
inverse scan variance as determined from the scan end data 
after filtering. A section of the mo survey map is shown in the 
top panel of Figure[3] Scan variances are coadded into a 'vari- 
ance map', which produces an estimate of the pixel variance 
over the survey area. 

4.2.2. Secondary Filtering 

Small-scale noise between rows of map pixels (visible in 
the top panel of Figure [3]l is further reduced by a destriping 
algorithm as follows. After the mo map is constructed, com- 
pact sources are located using the source extraction code de- 
scribed in the Source Paper, which is based on the SExtractor 
routine ( Berlin & Arnouts,1996^ . A second map (mi) is gener- 
ated identically to mo, with the exception that sources located 
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Fig. 3. — Maps illustrating the destriping process on QUaD 150 GHz data. 
All color scales are in MJy/sr, with saturated pixels within the survey region 
shown in white. Top: DC + slope filtered (mo) map. Bottom: Destriped map 

m the scan ends are masked during filtering. This process pre- 
vents discrete sources lying far from the bulk of the diffuse 
emission from influencing the inital polynomial filtering. 

The simple removal of a DC-level and slope from each 
scan produces a map which still exhibits striping due to at- 
mospheric 1 // noise. To suppress this noise, we construct a 
template for the sky signal, which is simply a smoothed ver- 
sion of the mi map. This template is then subtracted from the 
raw timestream, leaving data which is dominated by atmo- 
spheric noise. A 6th-order polynomial is then fit to the signal- 
subtracted timestream, and then subtracted from the original 
data. The timestream still contains the galactic signal of inter- 
est, but with the 1 // noise much suppressed compared to the 
simple DC-i-slope filtering described above. 

Maps made with the secondary-filtered data are termed the 
m2 maps — these exhibit improved noise properties without 
a heavy penalty in signal loss, as may be seen in Figure |3] A 
detailed discussion of the algorithm and its implementation is 
presented in Appendix [A| 

4.3. Absolute Calibration 

Absolute calibration is applied using conversion factors 
from [Brown etaL] ( [2009) . The maps used in the QUaD CMB 
analysis were cross-calibrated with the Boomerang experi- 
ment ( Masi et al.|200 6) to produce factors which convert from 
detector units of volts to thermodynamic units of /iK. These 
are estimated to be uncertain at the 3.5% level. The galac- 
tic maps are calibrated using the same conversion factors, and 
then rescaled to brightness units using 

dI = (dB/dT)2j3xdT, (1) 

where B is the Planck function, and dT and dl repre- 
sent thermodynamic and brightness fluctuations respectively. 
Througout the analysis in this paper, the central frequencies 
vo of the QUaD bands are loosely referred to as 100 and 
150 GHz. Assuming a spectrally flat source, numerical in- 
tegration of the QUaD bandpass presented in the Instrument 
Paper results in central frequencies of 94.5 and 149.6 GHz. 
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4.4. Final Maps 

Calibrated, destriped survey maps of Stokes /, Q, and U 
are shown in Appendix [C] Figures 19 21 in celestial coor- 
dinates (the native coordinate system for map processing) . 
The 100 GHz maps are hereafter referred to as /loo, Qioo, and 
f/ioo, and likewise at 150 GHz. All polarization maps follow 
the lAU convention, in which +Q is parallel to N-S and +U 
parallel to NE-SW (Ha maker & Bregman|1996| l. Converting 
to galactic coordinates (galactic longitude / and latitude b) re- 
sults in the maps shown in Figure [5] Total intensity maps are 
transformed from the native celestial coordinates to galactic I 
and b by linear interpolation of the map pixel values between 
the respective coordinate grids calculated with standard astro- 
nomical software packages. For polarization, we compute the 
angles between unit (RA,dec) vectors in the {l,b) basis at each 
point in the map — these angles are then used to project the 
polarized intensity onto the {l,b) basis. 
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Fig. 4. — Plots showing median sensitivity of the survey as a function of / 
and b for the fourth quadrant region data; plots for the third quadrant region 
are very similar. In all panels, black represents 100 GHz /, blue 150 GHz /, 
green 100 GHz Q, and red 150 GHz Q. The sensitivity to U is comparable 
to Q. Top row is for the full resolution (pixel size 0.02°) maps, while bottom 
row is for the coarse resolution (0.5° pixels) maps used for the analysis of 
diffuse structure in Section|5] 

The median sensitivity of the survey as a function of I and b 
is plotted in Figure|4]for the fourth quadrant region, where the 
rms values are determined from the variance maps; the sensi- 
tivity is almost identical for the third quadrant. The QUaD 
galactic plane survey is fairly uniform over ^ 50° in / and 
\b\ < 4° in each region, allowing detection and characteri- 
zation of diffuse and localized emission over ^ 800 square 
degrees of the low-latitude galaxy. 

At the native map resolution of 0.02°, much of the dif- 
fuse emission and hundreds of compact sources are detected 
to high significance; particularly bright regions such as the 
galactic center reach signal-to-noise ratios S/N > 500 in total 
intensity, with the bulk of the diffuse emission (within ^ 3° of 
the galactic plane) detected with S/N > 10 per map pixel. In 
Figure [4] we also show the sensitivity for 'coarse' resolution 
maps with pixel size of 0.5°; in these maps, the equivalent 
brightness sensitivity is a factor ^ 25 higher than the full res- 
olution maps. In polarization, there is significant diffuse po- 
larized emission in the 150 GHz fourth quadrant data at this 
resolution, while polarized signal at 100 GHz becomes appar- 
ent when degrading to 0.5° pixels. Note that when construct- 
ing coarse resolution maps, all timestream operations such 



as field-differencing, filtering and destriping are performed 
exactly as above. Locating point sources, which forms part 
of the destriping process, is done on the native 0.02° resolu- 
tion niQ and mi maps as above — the only difference between 
making the coarse and full resolution maps is the map pixel 
size which the processed timestream is binned into when con- 
structing the final m2 maps. 

The fainter diffuse 100 GHz signal implies that we are ob- 
serving emission with a positive spectral index a (I (x ly"), 
likely dominated by polarized thermal dust. As reported in the 
Source Paper, a small number of discrete sources are detected 
in the polarization maps, along with a polarized arc near the 
galactic center, and an extended polarized cloud. 

In addition to the signal maps, we also generate 'jackknife' 
maps in an idential manner to P09. The timestream data is 
split evenly in two ways, scan direction and time (i.e. first 
half of data against second half). For each jackknife, maps 
are constructed from each split of the data exactly as for the 
signal maps; the difference of pairs of maps from each split 
is then taken. These jackknife maps provide useful tests of 
possible systematics in the data or mapmaking process, and 
are generated for /, Q and U maps at both frequencies. 

5. PROPERTIES OF DIFFUSE EMISSION 

The QUaD data can be used to place constraints on inter- 
esting properties of the bulk galactic emission, namely the to- 
tal intensity spectral index, the polarization fraction, and the 
angle of polarization. Due to the low signal-to-noise of po- 
larization, we use the coarse resolution maps described above 
for the analysis; for convenience, the maps are interpolated 
to galactic coordinates as described in Section 4.4 On the 



coarse pixel grid, the effects of differing beam sizes between 
the two frequency bands become negligible compared to the 
effects of pixelization. 

Figure |5] shows polarization vectors from the coarse maps 
overlaid on full resolution total intensity maps for the \b\ < 
2.5° region of the survey. Pixels with polarized intensity 
signal-to-noise P/ap > 5 have been excluded, with the po- 
larized intensity 



and its error given by 



■ (7?) + (7?,+ 2(7? 



(2) 



(3) 



where Oq and afj are the variances in Q and U, and aqu is the 
covariance between Q and U map pixels. The 150 GHz data 
in particular show a high degree of alignment between the 
polarization (pseudo) vectors. In the statistical analysis that 
follows, pixels from the third and fourth quadrant maps are 
combined. Throughout, we ignore the contribution of primor- 
dial CMB anisotropics; the bulk of the analysis is performed 
within -1° < b < 1°, where the galactic emission is expected 
to be overwhelmingly dominant in both total and polarized in- 
tensity. In addition, pixels whose polarized emission is dom- 
inated by discrete sources rather than the diffuse background 
are also excluded. This flagging only includes the galactic 
center and a polarized cloud at I ^ 345°, b ~ 1.75°. To de- 
termine bulk emission properties, the analysis which follows 
makes use of the formalism presented in |Weiner et al. ( 2006 1, 
who define a generalised statistic L which accounts for er- 
rors in both coordinates e^. and e,,, and intrinsic scatter cr,. in a 
linear model y=A + Bx: 

,22,2, 2^ , (yi-(A + Bxdf 



(4) 
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Fig. 5. — Total intensity and polarization vector maps for third quadrant 100 GHz (first) and 150 GHz (second), and fourth quadrant 100 GHz (third) and 
150 GHz (fourth). The total intensity maps are full-resolution maps smoothed to the beam scale; the color scale is in MJy/sr. The polarization vectors are derived 
from coarse resolution maps with 0.5° pixels. Only polarization vectors for which P/ap > 5 are shown. 
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L is minimized using standard routines that return errors on 
the fit parameters A, B and tr,,; these errors are also verified 
by using bootstrap realizations of the data. Pixel values are 
the X, and y,, with their errors e^ i and e,./ determined from 
the variance maps. The intrinsic scatter term accounts for the 
variation in sky signal above that expected from the instru- 
mental noise alone. For the purposes of quick reference, a 
summary of the diffuse analysis is presented in Table [T] 



5.1. Total Intensity Spectral Index 

The spectral index in total intensity is calculated by min- 
imizing Equation |4] to find the slope B (with the intercept A 
held fixed) between the = /i5o_, and x, = /joo,,- pixels. Pix- 
els used in this analysis are restricted to those with \b\ < 1°, 
where the effects of data processing are smallest — see Ap- 
pendix |A3] The slope is simply related to the spectral index 
as 

ln(150/100)' ^' 

where the spectral index is calculated at the nominal QUaD 
center frequencies. These are assumed fixed, regardless of the 
spectral index of the underlying emission mechanisms — a 
more rigorous analysis would involve integrating each source 
model over the bandpass and re-calculating the central band 
frequencies, as the source spectral index can cause j/q to shift. 
We determined that for reasonable values of the source spec- 
tral index, fo varies by only a few percent, changing our re- 
sults insignificantly compared to absolute calibration errors. 
Converting to spectral index via Equation |5] the best fit pa- 
rameters are a, = 2.38 ±0.01 ±0.02, and a^i = 0.158 ±0.004. 
The first error on a/ is statistical, while the second is sys- 
tematic as estimated from signal-only simulations (see Ap- 



<a,>=2J5±0.0l; c =0.158+0.004 



pendix B.l 



To test for spectral index variations as a function of galactic 
latitude, the pixels are subdivided into rows of constant b, and 
the analysis above repeated. Figure |6] shows the results: The 
spectral index appears to flatten between +1° and -1°, moving 
from aj ^2.4 to aj ^ 1.8. In principle, the field-differencing, 
filtering and destriping processes could cause this measure- 
ment to be biased due to systematic effects. However, signal- 
only simulations described in Appendix B.l show that these 

7c and a scatter 
Variation of aj 



processes introduce a systematic shift of < 
of ^ 7% to the spectral index measurement 
with b in the QUaD data is therefore hkely a real property of 
the galaxy. 

The QUaD value of the spectral index between 100 and 
150 GHz is lower than that expected for dust alone (aj — 4), 
implying one (or more) additional emission component is 
present. Possible candidates for the extra component are syn- 
chrotron, free-free or spectral line emission. Synchrotron has 
a steeply falling spectral index (a., ^ -0.7), but is not expected 
to dominate near the plane on account of its large scale height. 
Alternatively, free-free emission has a flatter spectral index 
(ajf ^ -0.1) and is concentrated towards b = due to the 
collisional nature of the process. Given the expected spectral 
indices of these components, it is likely an excess of emis- 
sion in the 100 GHz band (rather than a 150 GHz deficit) that 
causes the relatively flat QUaD spectral index. Spectral line 
emission is a further possibility, with several authors report- 
ing large line contributions to the bolometric flux, though at 
higher frequencies (e.g. Groesbeck 1995 Comito et al. 2005 
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Fig. 6.-- Left: scatter plot of 150 GHz vs 100 GHz / pixels. The solid 
black line is the best fit mean slope B = /150//10O: with the broken black lines 
indicating the it 1 erg regions described by the intrinsic scatter of the data. B 
and (Tg are converted to spectral index a/ and intrinsic scatter CTc,, , with the 
best fit values and errors for these parameters shown above the plot. Right: 
Q/ as a function of galactic latitude b. The black dots show the mean value 
and its error, while the gray error bars indicate the intrinsic scatter (Tq, . 



[WyrowskTe t al. 2006). The two QUaD bands alone are insuf- 
ficient to separate these components; a discussion of spectral 
fitting in conjunction with the WMAP data is deferred to Sec- 
tion |5.3| while Section|5]4]discusses the relative contributions 
of synchrotron, dust and free-free as predicted from models. 

5.2. Polarization 

Two quantites are of interest from the polarization data: the 
polarization fraction, and the angle of polarization. The anal- 
ysis of polarization fraction proceeds similarly to the spec- 
tral index. This time, we search for the gradient and intrinsic 
scatter in a plot of Q or f/ against / via minimization of Equa- 
tion [4j the intercept is fixed at zero as before. Appendix |B.2| 
demonstrates that field-differencing, filtering and mapmaking 
processes bias the recovery of Q/I or U/I by ^0.1%. 

For the polarization fraction analysis, the jackknife maps 
are also used to test for contamination; the analysis is identi- 
cal to the signal data. Note that the polarization fraction for 
jackknife pixels is defined as 



_ Qjack 
jack ^signal 



(6) 



and likewise for U, where Isignai is the un-jackknifed total in- 
tensity map. Ideally there is no sky signal in the jackknife 
maps, and hence Equation |6]represents the polarization frac- 
tion that would have been observed if the sky had a true po- 
larization fraction of zero. 

The average 100 GHz polarization fractions are Q/I = 
1 .38 ± 0.06% with ct^// = 0.54 ± 0.02% for the intrinsic vari- 
ance, and U/I = -0.12 ± 0.05% with au/i = 0.51 ± 0.02%. 
At 150 GHz, 2/7= 1.68 ±0.04% with ctq// = 1.37 ±0.04%, 
and U/I = 0.27 ± 0.04% with a^/, = 1.22 ± 0.04%. Com- 
bining these results and subtracting the noise bias in P (i.e. 

Pdebias 

= \/P^-al), we find average polarization fractions of 
P/I= 1.38±0.08±0.1% at 100 GHz andP//= 1.70±0.06± 
0.1% at 100 and 150 GHz respectively, where the first error 
is random and the second due to systematic effects as deter- 
mined from signal-only simulations (Appendix |B.2| i. The cor- 
responding intrinsic scatter is simply computed as the quadra- 
ture sum of that from Q/I and U/I, and is 0.74 ±0.03 and 
1.83 ±0.06 at 100 and 150 GHz respectively. 

Figure It] displays the average Q/I and U/I polarization 
fraction or diffuse emission as a function of galactic latitude 
for the QUaD signal and jackknife data. The signal data 
shows polarization fractions which are close to constant with 
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TABLE 1 

Average Diffuse Emission Properties 



Property 1 


100 GHz 




150 GHz 






X 




X 


0".v 





2.35 ±0.01 ±0.02 


0.158 ±0.004 






P// (%) 


1.38±0.08±0.1 


0.74 ±0.03 


1.70 ±0.06 ±0.1 


1.83±0.06 


4> (deg) 


4.1± 1.3±5 


33.7± 1.0 


7.0±1.1±2 


29.7 ±0.9 



Note. — Summary of average diffuse emission properties from QUaD data. The symbols x and 
represent the mean and intrinsic scatter of each quantity in the leftmost column; where quoted, 
the first eiTor is statistical and the second systematic. 
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Fig. 7. — QUaD measurement of galactic-longitude averaged polarization 
fraction as a function of galactic latitude. Top row is 100 GHz Q/I (left) and 
U /I (right); the QUaD signal data are black points, with the scan direction 
and time jackknife are squares and crosses respectively (the jackknife data 
are offset in h from the signal data for clarity). Bottom row is for 150 GHz. 

b. The mean value in each jackknife b bin is expected to be 
consistent with zero for both Q/I and U /I, with the variance 
due to pixel noise only — Figure [7] shows that this is largely 
the case, with the jackknife data in each bin consistent with 
zero at the 2<t level or better 

The P// values of polariz ation fraction are so mewhat lower 
than the Archaeops result ( Benoit et al.||2004| l, who found a 
4-5% polarization fraction for \b\ <2" over the galactic longi- 
tude range 297 to 85° at 350 G Hz. Conversely, using WMAP 
3 year data, Kogut et al. (12007 ) found a 94 GHz polarization 



fraction closer to '--^ 1% averaged over a region including the 
QUaD survey between galactic latitudes -10 < b < 10°, ris- 
ing to 3.6% outside the P06 mask used in their analysis. A 
direct comparison of QUaD to the WMAP data is presented 
in Section |53] 

If the large-scale magnetic field is largely aligned in the 
plane of the galaxy, in the galactic coordinate system this 
translates to polarized emission predominantly in +Q. Fig- 
ure |7] demonstrates that this is observed by QUaD, though 
the U/I data at 150 GHz shows detected signal. This ob- 
servation may be alternatively quantified by directly comput- 
ing the polarization angle (p = 0.5tan"' iP 10) ^^'^ error 
d0 = Q.5{l+U^/Q^)-^A{U/Q), where 



d(f//e) 




(7) 



The effective function in Equation |4] is then minimized 
— this time, since we wish to know the mean and intrinsic 
variance of the distribution, the gradient term B is fixed 
at zero. The analysis is performed on the same pixels as 
polarization fraction, with the mean and intrinsic variance 



calculated for all pixels and in rows of constant b. Fig- 
ure |8] shows the results. The mean polarization angles are 
0100 = 4. 1 ± 1 .3° ± 5° and 0150 = 7.0 ± 1 . 1 ° ± 2°, with intrin- 
sic variance cr^, loo = 33. 7± 1.0° and cr^^igo = 29.7±0.9°. For 
0, the second error is the estimated systematic uncertainty in 
the pol ariza tion angle using signal-only simulations (see Ap- 
pendix B.2i. In bins of constant galactic latitude, and its 
intrinsic variance do not vary significantly. A weighted prob- 
ability distribution of is also shown at each frequency for the 
signal data, and all jackknife data combined. At both frequen- 
cies, a distinct peak is observed in the distribution close to the 
mean values calculated above, while the jackknifes are con- 
sistent with random numbers distributed uniformly between 
-90° < < 90°, as expected in the presence of no signal. 
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Fig. 8. — QUaD polarization angle; top and bottom are 100 and 
150 GHz respectively. Left: </> vs h in 0.5° bins of galactic latitude. Black 
points show the mean <)> in each bin, while black and gray error bars indicate 
the statistical error on the mean and the intrinsic scatter respectively. Right: 
weighted probability distribution Pr(4>) of the pixel values as a function of (f>. 
The black solid line corresponds to the data, while the broken gray line shows 
Pr((/)) for the combined jackknife maps. The QUaD data clearly show a peak 
in Pr(0) compared to the jackknife data, indicating the presence of coherent 
polarized galactic signal. 

The QUaD data indicate that while the galactic magnetic 
field is preferentially aligned parallel to the plane, there is 
significant additional scatter present. However, signal-only 
simulations (Appendix B.2 1 show that a substantial contribu- 
tion to the scatter in the angle may be present due to filtering 
and processing and effects. At 100 GHz this systematic scat- 
ter is comparable to the observed scatter, indicating we cannot 
reliably constrain cr^ at this frequency. 

5.3. Comparison to WMAP 

The closest comparison to the QUaD 100 GHz data is from 
WMAP W band, centi-ed on 9 4 GHz. For the follow ing anal- 
ysis, the 5-year WMAP data ( Hinshaw et al.||2009 i are used 
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<Q/I>=L13±0.32%;<U/I>=-0.24±0.27% 
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Fig. 9. — Left: QUaD 100 GHz and WMAP W bandpass spectral response in black and gray respectively. The QUaD band is substantially wider, with 
Af = 27 GHz compared to Au = 20 GHz for WMAP W band. Center: Example data points and spectral fit for a representative pixel in the WMAP and QUaD 
maps. WMAP data are gray crosses, QUaD data are black dots, and the broken black line indicates the best-fit two-component spectrum. The inset shows the 
distribution of spectral indices taken over map pixels used in the analysis. Right: Average polarization fraction as a function of galactic latitude. Gray is WMAP, 
black is QUaD 100 GHz for comparison; circles are Q, crosses are U. The points have been offset in b for clarity. 



to generate simulated timestream in the regions observed by 
QUaD. The timestream is processed in the same manner as the 
signal-only simulations described in Appendix |A] in which 
the field-differencing and filtering operations are performed 
exactly as for the QUaD data. 

The combination of QUaD and WMAP total intensity data 
allows constraints on the spectral index of emissive compo- 
nents in the galaxy. Using the QUaD and WMAP beam and 
pixelization functions, the WMAP Ka, Q, V and W band 
maps, and the QUaD survey, are convolved to the WMAP 
K-band resolution of 0.93°, and binned into the 0.5° pixels 
used in the coarse resolution QUaD maps. Pixel noise in the 
WMAP maps are determined from regions well away from 
the galactic plane, but are much smaller than the absolute cal- 
ibration uncertainties of the smoothed maps. Using the same 
pixel from each map, the data from each of the seven bands 
are fit to a two-component continuum model, the sum of two 
power laws in frequency i^: 



I{v)=A,u'''+AdV° 



(8) 



In this expression, a,, and are the synchrotron and dust 
spectral indices respectively. Note that the first component is 
only loosely termed as due to synchrotron; as discussed previ- 
ously, free-free is likely the second most dominant emission 
mechanism after dust at 60-100 GHz. To fit the data, the 
model is convolved across each bandpass to yield the average 
intensity in that band: 



/ I{v)T(v)dh 
J T{v)dv 



(9) 



where T{v) is the bandpass response as a function of fre- 
quency, as shown for WMAP W band and QUaD 100 GHz 
in the left panel of Figure [9] The against the data is cal- 
culated and minimized to find the best-fit parameters for the 
model in Equation |8] with Equation |9] evaluated for every set 
of trial parameters used in the minimization. 

The QUaD data is 25% brighter than WMAP at similar 
frequencies, indicating a discrepency between the two exper- 
iments which is further discussed in Section 15.41 The center 
panel of Figure |9] shows the fit spectrum for a single repre- 
sentative pixel in the data, and a scatter plot of a, against 
ad for the spectral fit to each map pixel. Taken over all pix- 
els, the average spectral indices are a., = -0.32 ± 0.03 and 



ctd = 2.84 ±0.03. The former is flatter than might be ex- 
pected for pure synchrotron (a ^ -0.7), indicating the pres- 
ence of free-free and/or dust, while the latter is lower than 
that predicted by the EDS dust model 8 {ctd.FDS = 3.5; see Ap- 
pendix [BT| . It is clear that the simple two-component model 
used is inadequate to describe the data; this statement remains 
true if the QUaD 100 GHz data is excluded. In fact, |Goldetal.| 
(2009 ) find that inside the plane (interior to the WMAP 5 year 
KQ95 mask), a ten-parameter model is insufficient to fully de- 
scribe their data, and therefore it may be optimistic to expect 
simple models to account for all emission along lines of sight 
close to the plane. 

The right panel of Figure |9] compares the polarization 
fraction, averaged along galactic longitude, between the un- 
smoothed WMAP W band maps, and QUaD 100 GHz— both 
maps are binned at 0.5° resolution. The plot demonstrates that 
QUaD is consistent with the WMAP observations to within 
1(7, in both Q and U, with WMAP mean values of Q/I = 
1.13±0.32% and t/// = -0.24 ± 0.27% (P/I= 1.1±0.4%), 
compared to 1.38 ±0.06% and -0.12 ±0.05% for QUaD 
(P/I = 1.38 ±0.08%); agreement between the two datasets 
is also found as a function of galactic latitude. 

5.4. Comparison to Emission Models 

Synchrotron and dust models are provided by |Finkbeiner| 
( |2001| ) and [Finkbeiner et al.| \\999 \ (we use dust model 8 
in the latter case); these can be extrapolated into the QUaD 
bands and compared to the observations. At frequencies 
^100 GHz, free-free emission is also expected to contribute. 
The Ha sky template provided by |Finkbe iner (2003) is used to 
generate a free-free emission map at the QUaD frequencies; 
we follow jSchafer et al.| ( ,2006) , who convert the Ha tem- 
plate from Rayleigh units into /LtK using the formula provided 
by |Valls-Gabaud| ( ,199"8l l: 



14.0 



0.317 



290 r V \~ 



(10) 

where Ana are the template pixel values in Rayleighs, Tp is 
the plasma temperature (assumed to be lO'^K), v is the central 
observing frequency of the QUaD bands, an d gff is the free- 
free Gaunt factor as calculated in Finkbeiner| ( |2003| ). The tem- 
perature maps are converted to brightness units as in Equa- 
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tion[Tl 

Signal-only simulations of synchrotron, dust and free-free 
are used to generate field-differenced and filtered maps of 
these components, which are then summed to produce a 
model sky at QUaD frequencies. The limiting factor in resolu- 
tion is the synchrotron model, which is derived from 405 MHz 



maps of [Haslam etaL] (1981] |1982[ ) with a beam FWHM of 
~ 1°. QUaD and model maps are binned into 0.5° pixels and 
smoothed to 1° resolution to match the synchrotron model. 

In Figure 10 we plot the pixel values between QUaD and 
the model predictions against each other (left panels), and 
the mean ratio of pixel values between QUaD and the mod- 
els R = iQUao/Imodei in bins of b (right panels). R is the 
gradient in a plot of IquoD against I„,„der, this is calculated 
for all pixels and as a function of b by minimizing Equa- 
tion]?] holding the intercept fixed as in Section [5T| with in- 
trinsic scatter fitted simultaneously. The models are split 
into three compositions: dust alone, dust-nsynchrotron, and 
dust-HsynchrotronH-free-free, with comparisons to QUaD data 
made for each. 
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Fig. 10. — Left panels: QUaD pixel values plotted against signal-only 
simulations of the model including dust, synchrotron, and free-free as de- 
scribed in the text; top is 100 GHz, bottom is 150 GHz. Right panels: Ratio 
of QUaD map pixel values to model predictions at 100 (top) and 150 GHz 
(bottom), as a function of galactic latitude. Square symbols compare QUaD 
to the FDS dust model, crosses are dust+synchrotron, and dots correspond to 
dust+synchrotron+free-free. 



At 100 GHz, QUaD is brighter than the FDS dust-only pre- 
diction by a factor 1 .66 ± 0.01, decreasing to 1 .49 ± 0.01 and 
1.48 ±0.01 as synchrotron and free-free models are added; 
the intrinsic scatter is ^ 0.1 for all models. This extra signal 
at 100 GHz is hereafter referr ed to as the 'QUaD excess' — 
we note that |Gold et al.| ( |2009] l also find an excess of observed 
W band emission over that predicted by the Finkbeiner et al. 
(1999) models over most of the sky. At 150 GHz the model 
predictions are in better agreement with QUaD; R ranges from 
0.962 ± 0.007 (dust only) to 0.942 ± 0.006 (all model com- 
ponents), with the intrinsic scatter ~ 22%. The ratio R only 



changes by 1-2% comparing QUaD to the model with dust 
only or all components at this frequency; this is to be ex- 
pected since dust dominates the emission at 150 GHz. At 
both frequencies, R does not vary strongly with galactic lat- 
itude, though care should be taken when interpreting these 
results since the maps are smoothed to 1° resolution, and thus 
the data points are highly correlated between -1° < b < 1°. 

It should be noted a factor ~ 2 uncertainty exists in the 
conversion factor from Rayleigh units to antenna tempera- 
tures which we apply to the Finkbeiner (2003) free-free maps, 
which could contribute to the QUaD 100 GHz excess. How- 
ever, Figure[TO]shows that the addition of a free-free template 
to a dustH-synchrotron model changes the ratio of QUaD to 
model pixel values by less than 1%; thus even a factor 2 un- 
derestimation of the free-free model calibration is insufficient 
to account for the excess emission. 

Another possible explanation of the QUaD excess is molec- 
ular line emission: large fractional contributions from line 
emission have been meas ured at higher frequencies towards 
star-forming regions (e.g. Groesbeck||1995 [Nummelin et al. 
1998 1; these can range from 10-65% of the bolometric in- 
tensity. On the other hand, the COBE FIRAS instrument has 
detected much lower line emission contributions (<1%) over 



large patches of sky (e.g ., [Wright et al.||1991] [Bennett et al. 
TWi Tixsen et al.|1999| l. 

If the line emission interpretation is correct, one might ex- 
pect the QUaD 100 GHz data to be brighter than WMAP 
W band on account of the wider bandwidth admitting more 
lines (see Figure |9]l, and the WMAP data to be brighter than 
the models which do not include line emission at all. This 
trend is indeed observed: QUaD 100 GHz data is a factor 
1 .48 brighter than the model with dust, synchrotron and free- 
free included, and QUaD 100 GHz is also ~ 25% brighter 
than WMAP W band. WMAP is therefore some 23% brighter 
than the combined models for continuum emission, providing 
independent evidence that these models are insufficient to de- 
scribe the data close to the plane of the galaxy. 

5.5. Tests of the Spectral Line Hypothesis 

Given the variation seen in the literature on spectral line 
contributions at higher frequencies, the hypothesis that line 
emission causes the QUaD excess is subjected to the follow- 
ing additional tests. 

5.5.1. Spectral Comparison Using FIRAS 

The FIRAS instrument aboard the COBE satellite provides 
absolute spectral measurements covering the QUaD 100 GHz 
band, at a resolution of ^ 7°. A direct comparison between 
QUaD and FIRAS is not possible since the FIRAS beam 
width is half the maximum width (in R.A.) of the QUaD 
survey, and hence the QUaD maps cannot be smoothed to 
the same spatial resolution. Neither can spectral discrimina- 
tion be used to compare the FIRAS data as filtered through 
the QUaD and WMAP bandpasses, since even at maximum 
spectral resolution of 3.4 GHz, the frequency sampling is too 
sparse to resolve the difference between QUaD and WMAP 
where there is no spectral overlap. However, the FIRAS data 
can be used to obtain an upper limit on the emission expected 
in the QUaD 100 GHz band, providing a useful consistency 
check of the QUaD absolute calibration. We restrict ourselves 
to pixels less that 5° from the galactic plane, and subtract the 
best-fit CMB monopole blackbody spectrum from each FI- 
RAS frequency channel. Other filtering effects such as poly- 
nomial subtraction are not included because the large FIRAS 
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beam smooths the galactic signal to greater galactic latitudes 
than covered in the QUaD survey. Polynomial subtraction of 
the signal lying at the QUaD scan ends would then reduce the 
amplitude below the level of signal loss due to filtering in the 
QUaD survey itself. 

The signal-to-noise at ^100 GHz is low in the FIRAS 
data, so an upper limit is obtained by finding the 95th per- 
centile of the pixel distribution for each frequency channel. 
This 'spectrum' of 95th percentiles is then integrated over the 
normalized QUaD 100 GHz bandpass, resulting in an upper 
limit of 5.3 MJy/sr at 95% confidence. This limit is consis- 
tent with the QUaD data; Figure |5] shows that at 100 GHz 
the peak brightness is ^ 7MJy/sr near the galactic center, 
which would be reduced by the ratio of beam areas (a fac- 
tor ^ (7 X 60/5)^ = 7056) once beam smoothing is taken into 
account. Therefore the FIRAS data cannot isolate the QUaD 
excess as being due to spectral lines or an absolute calibration 
mismatch. 

5.5.2. Spectral Comparison Using CO 1-0 Transition Maps 

The most prominent line near the lower QUaD band is the 
1-0 CO transition at ^ 115 GHz. This has been mapped 
over the inner galactic plane by Da me et ST] ( |2001[ ) and refer- 
ences therein; the survey has spatial resolution 6', very close 
to QUaD at 100 GHz. As seen in Figure |9] the QUaD band 
response at 1 15 GHz is small (a factor of ^ 1000 smaller than 
the peak in fact). The CO maps are restricted to the QUaD sur- 
vey boundaries, and the antenna temperature units converted 
to MJy/sr, first by converting antenna temperatures to ther- 
mal temperature units for each frequency channel, and then 
by averaging the CO data in frequency over the spectral band- 
width, and finally averaging over the QUaD bandpass. We 
find a peak CO contribution of 1 0~^ MJy / sr, or approximately 
2% of the typical brightness of a WMAP pixel in the QUaD 
survey region, and is therefore insufficient to account for the 
^ 25% excess seen by QUaD over WMAP. 

5.5.3. Spectral Line Check Using Polarization 

Spectral lines are not expected to emit polarized radiation, 
so any contribution can be tested by comparing the unpolar- 
ized and polarized data. 



The left panel of Figure 1 1 shows a plot of QUaD 150 vs. 
100 GHz data in / and Q. The slope in / is clear, and is appar- 
ently traced by Q\ this indicates that line emission does not 
contribute significantly to the lower QUaD band, as otherwise 
a steeper slope would be expected in polarization. However, 
the signal-to-noise at 100 GHz is low: Also plotted in Fig- 
ure 11 are the pixel values from the scan-direction jackknife, 
which demonstrate that instrumental noise contributes signif- 
icantly to the 100 GHz data, and therefore statistical uncer- 
tainties could bias the polarized spectral index measurement 
towards lower values. 

Greater signal-to-noise can be acheived by using the galac- 
tic longitude-averaged data. Neither dust nor synchrotron are 
expected to vary polarization fraction with frequency, yet as 
seen in Table[T]and Figure 1 1 the observed polarization frac- 
tions differ between QUaD bands. Taking the average values 
of P/I at each frequency, the unpolarized spectral index, and 
combining the measurement uncertainties and intrinsic scat- 
ter into a single error term, we find ap = 2.81;!;q25. This is 
discrepent with the total intensity spectral index at the ^ 2a 
level, and is thus an inconclusive test on the existence or not 
of additional emission components. 



As a final test, the QUaD 100 GHz maps are compared to 
WMAP W band in Stokes Q on a pixel-by-pixel basis. The 
right panel of Figure [TT] shows such a plot of this test; both 
data sets are noisy and a statistically significant measurements 
of the gradient is not possible. Instead, we simply take the 
mean of each set of pixels and calculate the ratio, finding 
{Q)QUaDl {Q)wMAP = 1.4 ±0.5. Such large errors render the 
data insufficient to measure an excess of QUaD polarization 
over WMAP, and thus cannot verify or falsify a contribution 
due to Une emission or an absolute calibration mismatch. 
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Fig. 11.— Left: Scatter plot of QUaD 150 vs. 100 GHz data. Black 
points are for unpolarized intensity /, red are for Q, and blue are for the 
scan-direction Q jackknife to indicate the level of noise present in the QUaD 
100 GHz Q data. The polarization data have been multiplied by 50 for visual 
display purposes only. Right: Scatter plot for QUaD 100 GHz data against 
WMAP W band for Stokes Q. The solid lines indicate the mean value of each 
data set. 



6. CONCLUSIONS 

We present the QUaD survey of the Milky Way in Stokes 
/, Q and U, with a resolution of 5 (3.5) arcmin at 100 
and 150 GHz respectively. The survey covers two regions, 
110° < RA < 190° and 210° < RA < 290°, both in the 
dec range -60.5° < dec < -26.5°, corresponding to approx- 
imately 245-295° and 315-5° in galactic longitude, and 
-4 < < 4° in galactic latitude — a total of ^ 800 square 
degrees. 

Degrading the map resolution to 0.5° pixels, the average 
spectral index of diffuse emission is a/ = 2.35 ± O.Ol(stat) ± 
0.02(sys), assuming nominal band centers of 94.5 and 149.6 
GHz. This value of spectral index is fla tter than that expected 
from dust alone — |Gold et al.| P009| l demonstrate that the 
WMAP 5 year data only weakly constrains the dust spectral 
index, but use a prior range 3.5 < < 5, indicative of the 
expectation for this component. The low QUaD-only value is 
interpreted as evidence for additional emission components in 
the lower frequency QUaD band. 

A direct comparison to WMAP 5 year W band data shows 
the QUaD 100 GHz maps are on average ^25% brighter Fit- 
ting a two-component continuum model to all pixels in the 
five WMAP and two QUaD bands results in constraints of 
a, = -0.32 ± 0.03 and = 2.84 ± 0.03. The first is attiibuted 
to a composite of synchrotron, free-free and dust expected at 
^ 100 GHz close to the galactic plane, with the second inter- 
preted as mostly dust. However, the fit is poor for such a sim- 
ple model and more emission components would be required 
to fully explain the data. 

Similarly to the spectral index determined from QUaD 
alone, is lower than the expectation from available models. 
A composite model of dust, synchrotron and free-free emis- 
sion underestimates the brightness at 100 GHz, where QUaD 
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observes on average a factor 1.48 ±0.01 more signal. One in- 
terpretation of this effect is molecular line emission. This pos- 
sibility is discussed further in Section [53] where a variety of 
cross-checks indicate that the 115 GHz CO 1-0 line is unlikely 
to be the main cause of the 100 GHz excess, but that the QUaD 
data are consistent with absolute spectral measurements from 
the COBE FIRAS instrument. At 150 GHz the agreement is 
better, with an average pixel ratio of 0.942 ± 0.006 between 
QUaD and the models. 

The QUaD data allow measurement of the polarization 
fraction in both bands — the results quoted here are taken 
from galactic coordinates maps using the lAU convention for 
Stokes parameters Q and U. Analysis in the Source Paper 
shows that few compact objects have measurable polariza- 
tion and thus we assume the dominant source of the polar- 
ized emission studied here is diffuse. On average, Q/I = 
1.38 ±0.06% and 1.68 ±0.04% at 100 and 150 GHz, with the 
equivalent averages for C/// being -0.12 ± 0.05% and 0.27 ± 
0.04%. The intrinsic scatter on these quantities were found 
to be typically - 0.5 1.3)% at 100 and 150 GHz respec- 
tively, reflecting fluctuations in polarization fractions at dif- 
ferent positions in the galactic plane. Signal-only simulations 
indicate that the systematic error on polarization fraction is 
of order 0.1% at both frequencies. Measurements of Q/I and 
U/I also possible as a function of galactic latitude b within 
\b\ < 1°, and show evidence of small deviations from the aver- 
age polarization fractions quoted above, but within the range 
allowed by the intrinsic scatter. Combining Q/I and U/I, we 
find total polarization fractions P/I = 1.38 ±0.08 ±0.1% at 
100 GHz and 1.70 ±0.06 ±0.1% at 150 GHz, where the first 
error is random and the second systematic. The intrinsic scat- 
ter at these frequencies is 0.74 ± 0.03 and 1 .83 ± 0.06. 

Comparing the QUaD polarization fraction to that from 
WMAP 5 year data, agreement is found between the two 
datasets at QUaD 100 GHz and WMAP W band, the lat- 
ter giving an average polarization fraction 1.1 ±0.4%. The 
polarization fraction measurements reported here provide en- 
couragement that large areas of sky may be useful for prob- 
ing inflationary cosmology with CMB polarization B-modes 
at frequencies above 100 GHz. 

The mean angle of polarization close to the plane measured 
by QUaD is is <</)>= 4. 1 ± 1 .3 ±5° at 100 GHz, with <(j>>= 
7.0±1.1±2° at 150 GHz, where the first error is random and 
the second systematic. Although the data indicate QUaD has 
detected intrinsic scatter in the distribution of (p, the amount 
of scatter is consistent with that introduced by filtering and 
map processing effects. The observations therefore provide 
evidence of the large-scale alignment of the galactic magnetic 
field. 

Extensive tests could not definitively isolate the cause of 
the QUaD excess signal observed over WMAP W band and 
continuum emission models at 100 GHz. The QUaD data is 
consistent with an unpolarized intensity upper limit derived 
from FIRAS data, ruling out a large QUaD absolute calibra- 



tion error A second test of the absolute calibration is calcu- 
lated from the ratio of average Stokes Q in QUaD and WMAP 
pixels; we find (Q) quoD / {Q)wmap = 1.4±0.5, an excess over 
/ on average, but statistically consistent with the mean total 
intensity excess of ^ 1.25. Low signal-to-noise in polariza- 
tion at 100 GHz prevented a direct determination of the polar- 
ized spectral index; a larger value than total intensity would 
be evidence for line emission in the 100 GHz band. Using a 
combination of the average polarization fraction and the un- 
polarized spectral index, we determined ap = 2.81;!;o within 
2(T of ai and thus inconclusive regarding molecular line emis- 
sion. Maps of the CO 1-0 transition ( |Dame et al.|2001| ) mul- 
tiplied through the QUaD bandpass placed an upper limit on 
the contribution of this molecular line of 10"^ MJy/sr, which 
at ^ 2% of the typical brightness of a WMAP W band pixel is 
insufficient to account for the QUaD excess. We conclude that 
higher signal-to-noise measurements of the polarized galactic 
emission at 100 GHz are required to resolve the QUaD ex- 
cess, an d should be provid ed in the near future by the Planck 
satelHte ( |Villa et al.|2002] l. 

This paper focussed on the properties of diffuse emission 
(typically 0.5° scales and larger); however, the QUaD galactic 
plane survey contains information on scales down to 5 and 3.5 
arcminutes at 100 and 150 GHz . The small-scale properties 
of emission via discrete sources are explored in the Source 
Paper, where a variety of objects have been detected, such as 
ultra-compact HII regions and supernova remnants. 

ACKNOWLEDGEMENTS 

This paper is dedicated to the memory of Andrew Lange, 
who gave wisdom and guidance to so many members of 
the astrophysics and cosmology community. His presence 
is sorely missed. We thank our colleagues on the Bl- 
CEP experiment and Dan Marrone for useful discussions. 
QUaD is funded by the National Science Foundation in the 
USA, through grants ANT-0338138, ANT-0338335 & ANT- 
0338238, by the Science and Technology Facilities Coun- 
cil (STFC) in the UK and by the Science Foundation Ire- 
land. The BOOMERanG collaboration kindly allowed the 
use of their CMB maps for our calibration purposes. MZ 
acknowledges support from a NASA Postdoctoral Fellow- 
ship. PGC acknowledges funding from the Portuguese FCT. 
SEC acknowledges support from a Stanford Terman Fellow- 
ship. JRH acknowledges the support of an NSF Graduate 
Research Fellowship, a Stanford Graduate Fellowship and a 
NASA Postdoctoral Fellowship. YM acknowledges support 
from a SUPA Prize studentship. CP acknowledges partial 
support from the Kavli Institute for Cosmological Physics 
through the grant NSF PHY-01 14422. EYW acknowledges 
receipt of an NDSEG fellowship. We acknowledge the use of 
the Legacy Archive for Microwave Background Data Anal- 
ysis (LAMBDA). Support for LAMBDA is provided by the 
NASA Office of Space Science. 



REFERENCES 



Ade, R et al. 2008, ApJ, 674, 22 
Bennett, C. L. et al. 1994, ApJ, 434. 587 
Benoit, A. et al. 2004, A&A. 424, 571 
Bertin, E., & Arnouts, S. 1996, A&AS, 1 17, 393 
Brown, M. L. et al. 2009, ApJ, 705, 978 

Comito, C, Schilke, R, Phillips, T. G., Lis, D. C, Motte, E, & Mehringer, D. 

2005, ApJS, 156, 127 
Dame, T. M., Hartmann, D., & Thaddeus, R 2001, ApJ, 547, 792 



Dunkley, J. et al. 2009, in American Institute of Physics Conference 
Series, Vol. 1141, American Institute of Physics Conference Series, 
ed. S. Dodelson. D. Baumann, A. Cooray, J. Dunkley, A. Fraisse, 
M. G. Jackson, A. Kogut, L. Krauss, M. Zaldarriaga, & K. Smith , 222-264 

Finkbeiner, D. P. 2001, private communication 

. 2003, ApJS, 146, 407 

Finkbeiner, D. R, Davis, M., & Schlegel, D. J. 1999, ApJ, 524, 867 
Fixsen, D. J., Bennett, C. L., & Mather, J. C. 1999, ApJ, 526, 207 
Gold, B. et al. 2009, ApJS. 180, 265 



QUaD Galactic Plane Survey 



13 



Input FDS model 



model + noise 



destriped timestream 



> 
S 





95 

azimuth (deg) 

Fig. 12. — Destriping process for one half-scan of data for the central 150 GHz pixel. The dashed graylines show the regions of the scan ends used for 
DC-level plus slope filtering. Left: Simulated input FDS model 8 signal. Center: Input model plus simulated noise before (black) and after (gray) initial filtering. 
Right: Destriped timestream — the recovered timestream is an improved estimate of the true sky signal compared to the initial filtering. 



Groesbeck, T. D. 1995, PhD thesis, AA(Califomia Institute Of Technology.) 

Hamaker, J. P., & Bregman, J. D. 1996, A&AS, 117, 161 

Haslam, C. G. T., Klein, U., Salter, C. J., Stoffel, H., Wilson, W. E., Cleaiy, 

M. N., Cooke, D. J., & Thomasson, P 1981, A&A, 100, 209 
Haslam, C. G. T., Salter, C. J., Stoffel, H., & Wilson, W. E. 1982, A&AS, 47, 

1 

Heiles, C. 1996, ApJ, 462, 316 

Hildebrand, R. H., Davidson, J. A., Dotson, J. L., Dowell, C. D., Novak, G., 

& Vaillancourt, J. E. 2000, PASP, 1 12, 1215 
Hinderks, J. et al. 2009, ApJ, 692, 1221 
Hinshaw, G. et al. 2009, ApJS, 180, 225 
Hu, W., & White, iVl. 1997, New Astronomy, 2, 323 

Jones, W. C, Bhatia, R., Bock, J. J., & Lange, A. E. 2003, in Presented at the 
Society of Photo-Optical Instrumentation Engineers (SPIE) Conference, 
Vol. 4855, Millimeter and Submillimeter Detectors for Astronomy. Edited 
by Phillips, Thomas G.; Zmuidzinas, Jonas. Proceedings of the SPIE, 
Volume 4855, pp. 227-238 (2003)., ed. T. G. Phillips & J. Zmuidzinas, 
227-238 

Kogut, A. et al. 2007, ApJ, 665, 355 

Lazarian, A. 2003, Journal of Quantitative Spectroscopy and Radiative 

Transfer, 79, 881 
Leitch, E. M. et al. 2002, ApJ, 568, 28 
Masi, S. et al. 2006, A&A, 458, 687 



Netterfield, C. B. et al. 2009, ArXiv e-prints, 0904.1207 

Nummelin, A., Bergman, P., Hjalmarson, A., Friberg, P., Irvine, W. M., 

Millar, T. J., Ohishi, M., & Saito, S. 1998, ApJS, 117, 427 
Olmi, L. et al. 2009, ArXiv e-prints, 0910.1097 
Ponthieu, N. et al. 2005, A&A, 444, 327 

Prunet, S., Sethi, S. K., Bouchet, F R., & Miville-Deschenes, M. 1998, A&A, 
339, 187 

Pryke, C. et al. 2009, ApJ, 692, 1247 

Schafer, B. M., Pfrommer, C, Hell, R. M., & Bartelmann, M. 2006, MNRAS, 

370, 1713 
SchuUer, F et al. 2009, A&A, 504, 415 

Valls-Gabaud, D. 1998, Publications of the Astronomical Society of 
AustraHa, 15, 111 

Villa, F. et al. 2002, in American Institute of Physics Conference Series, Vol. 

616, Experimental Cosmology at Millimetre Wavelengths, ed. M. de Petris 

& M. Gervasi, 224-228 
Weiner, B. J. et al. 2006, ApJ, 653, 1049 
Wright, E. L. et al. 1991, ApJ, 381, 200 

Wyrowski, F, Heyminck, S., Gusten, R., & Menten, K. M. 2006, A&A, 454, 
L95 

Zweibel, E. G., & Heiles, C. 1997, Nature, 385, 131 



APPENDIX 

A. FURTHER DETAILS ON DESTRIPING ALGORITHM 
A. 1 . Map Destriping 

When constructing the initial maps ;«o and mi, the choice of polynomial filtering order is a trade-off between increased at- 
mospheric noise reduction (higher order) and reduced filtering of the galactic signal of interest (lower order). A simple DC- 
level-i-slope filter function is used in the QUaD survey because we are interested in the emission properties on both small and 
large angular scales, both of which are suppressed or corrupted by a higher order filter function. However, this choice results in a 
larger atmospheric noise contribution than would be present had a higher order polynomial been used — the mo map in Figure |3] 
exhibits large row -to-row striping as a result, which we suppress as follows. After filtering, l/f noise is largely uncorrected 
between rows of pixels, while galactic structure is strongly correlated on these angular scales due to its intrinsic structure and 
beam smoothing. Smoothing the mi maps with a circularly symmetric gaussian kernel of width 4 x (Ji,ea„, at each frequency 
reduces striping between pixels; the smoothed map is treated as a template map of the sky, m,. From m,, 'signal only' timestream 
d, is interpolated and subtracted from the original data d: 

dn = d-d,. (Al) 

The 'signal-subtracted' timestream d„ is now dominated by atmospheric noise with little galactic signal present, and is fit with a 
higher order polynomial to measure the atmospheric modes (a 6th-order polynomial is used in the QUaD survey). 

Using the resulting polynomial coefficients p„, a polynomial t/^,, which represents atmospheric modes is subtracted from the 
original unfiltered timestream 

ds = d- dp^ , ( A2) 

and the filtered data is coadded into the map. Since dp^ largely measures atmospheric l/f modes, the noise in maps made from 
ds (denoted 1112 in the map ordering of our algorithm) is considerably whiter than the mo maps. The improvement may be seen in 
simulated timestream in Figure [12] and in the resulting maps in Figure [13] the procedure is repeated independently for /, Q and 
U maps at both frequencies. 

Some subtleties are present in this method; since m, is a smoothed version of the sky as observed by QUaD, subtracting it from 
the timestream introduces residuals near the locations of bright point sources (effectively, subtracting two gaussian functions of 
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equal cross-sectional area but differing widths). These residuals can be large, particularly near the galactic center, influencing the 
Pn, causing spurious filtering residuals in the 1112 maps. This effect is reduced by locating bright sources in mi as described in 
the Source Paper and then masking them as described below. Since 1 // noise in the mi total intensity maps can be spuriously 
detected as false sources, a high signal-to-noise threshold of S/N > 10 is used in the / map. For Q and U the noise is considerably 
whiter, so a threshold of S/N > 3 is used. 

With the locations of bright sources known, two methods are used to reduce their impact on higher-order filtering. First, the 
source pixels are replaced with a local median and the smoothed template map m, is computed — this reduces the amount of 
source power smeared out by the smoothing, lowering the corresponding residuals in d„. Second, the source locations are masked 
with a conservative radius of 6ai,eam when performing the 6th-order polynomial fit. 



After these steps, the data is filtered as in Equation A2 and coadded into the map, with post-filter inverse scan variances used as 
weights. Note that the variances for the destriped data are computed over the entire scan after template subtraction and filtering, 
and with point sources masked, i.e cr^ = var{d„-dpj. 

The mapping algorithm implemented in this paper can be summarized as follows: 

Step 1: Construct otq map from timestream filtered with DC-level+slope determined from scan ends. 
Step 2: Locate sources using the metliod described in the Source Paper, mask sources and and repeat Step 1 above to 
give mi. 

Step 3: Smootli mi with gaussian kernel of a, = 4x abeam, vvith bright source pixels replaced by local median - this is 
the template map mt. 

Step 4: Interpolate 'signal-only' timestream d, from m,, and subtract from the original data d to give an approximate 
'noise-only' timestream d„ . 

Step 5: Fit d„ with high order polynomial, masking sources located in Step 2. This step measures the 'noise-only' 
polynomial coefficients p„. 

Step 6: Filter the original timestream d using a polynomial with coefficients p„, and coadd the filtered half-scans into 
the map using post-filter inverse half-scan variances as weights; this produces the final m2 maps. 

Though the destriping algorithm is implemented entirely in timestream/image space, it has a direct interpretation as a linear 
filter in Fourier space as discussed in Section[A!4| 



A. 2. Test of Mapping Algorithm 

To test the algorithm described above, signal-only simulations of the FDS dust model 8 evaluated at the QUaD center fre- 
quencies are generated as mock sky signal. We wish to add realistic noise to the simulations, but since the galactic signal can 
dominate the timestream, we cannot directly take the power spectrum of the data and use it to regenerate noise as in P09 — doing 
so results in noise strongly correlated with regions of bright galactic emission. To mitigate this effect, the sky signal (as estimated 
from the destriped maps) is first removed from the data, and the resulting sky-subtracted timestream processed in the same way 
as P09: continuous segments of data are Fourier transformed, binned, and the covariance matrix of the fourier modes taken be- 
tween all channels in each bin separately. Noise timestream is then regenerated by mixing uncorrected random numbers with the 
Cholesky decomposition of this covariance matrix. Taking the inverse Fourier transform yields simulated noise timestream with 
the observed degree of covariance as the signal-subtracted data. P09 demonstrates that this process yields simulated noise which 
is indistinguishable from the real. The simulated noise is added to the signal-only timestream, coadded into 0.02° maps, and 



destriped according to Section A.l Figure [13] shows the signal-only input map m,, the DC-nslope filtered map mo, the destriped 
map m2, and the residual map m2-ms. Qualitatively, it is clear that the algorithm is effective at suppressing 1// noise, and 
introduces small residuals compared to the galactic signal of interest. 

A. 3. Effects of Map/Timestream Processing 

The mapmaking algorithm described above results in a loss of signal due to the main stages of processing: field-differencing 
and destriping. To test the effects of each, the same signal-only simulations of FDS model 8 as in Section [A!2] are used, comparing 
pixel values at each (cumulative) stage of processing of the input maps. 



The left panels of Figure 14 show the median fractional change in pixel values of mo maps, before and after field-differencing. 
At the mean 3(j noise level in the QUaD maps, field-differencing reduces the signal by 1-7% at 100 GHz, and 5-25% at 
150 GHz, depending on dec. The declination dependence arises from the fact that our 15° azimuth scans correspond to a smaller 
RA range at lower dec, scaling as cos(dec). Therefore the low dec trail field scans lie closer to the galactic plane, and so subtract 
out more sky signal when the lead and trail fields are differenced. The fractional loss in signal is smaller for brighter pixels close 
to the galactic plane, with less than 5% reduction above 0.35 (1.3) MJy/sr at 100 (150) GHz over all dec. 



A similar analysis is shown for the effect of destriping in the right hand panels of Figure 14 Since the destriping method only 



reduces uncorrected noise between rows of pixels, it is not dec-dependent. Destriping does introduce fluctuations in the pixel 



values, shown by the error bars in Figure 14 The fractional rms (relative to the signal) is at most 15% for pixels with amplitude 
above the 3a noise level at both frequencies, and decreases with increasing signal amplitude — the rms is < 5% above 0.5 (0.8) 
MJy/sr at 100 (150) GHz. 

A.4. Fourier Plane Interpretation of Map Destriping Algorithm 

Experiments such as QUaD which use fixed elevation scans suffer 1// noise predominantly in the scan direction. The 1// 
noise appears as a vertical band in the two-dimensional power spectrum Cu = /(u)/(u)*, where / is the Fourier transform of the 
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Fig. 13. — Maps illustrating the major stages of the mapmaking process. All colorscales are in MJy/sr, with saturated pixels within the survey region shown 
in white. Top left: Signal-only FDS model 8 prediction at 1 50 GHz. Top right: signal plus noise simulation of the same model after DC-level -I- slope filtering, 
i.e. mo stage map. Bottom left: Recovered mi inap after destriping. Bottom right: Difference between input and destriped maps. The color stretch has been 
reduced b y a fa ctor 5 to enhance the residuals; the noise consists mainly of white noise, and large-scale modes which are removed before the destriping in Step 4 
of Section lAJ] 
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Fig. 14. — Effect of field differencing and destriping on map pixel values, for FDS model 8 signal-only simulations. Top and bottom are 100 GHz and 
150 GHz respectively. Left: Mean fractional change in pixel values of mg map, before and after field differencing, as a function of input pixel value. Lines in 
different shades of gray correspond to different bins in dec, as denoted in the legend. Right: Mean fractional change in pixel values before and after destriping 
process (i.e. comparing m\ and mi maps). The error bars correspond to the standard deviation in each bin in logjo/ (MJy/sr). 



total intensity map, and u = (m, v) is the wavevector. The leftmost panels of Figure 15 show a signal-only map and Cu from an 
FDS simulation of a subsection of the survey; an example signal-nnoise mo map and corresponding Cu are displayed in the center 
panels, where the vertical noise band is clearly visible in the 2D power spectrum. Polynomial filtering over the entirety of each 
scan digs a 'trench' into the signal and noise in the vertical band: the width of the trench in Fourier space, Wy„ is determined by 
the polynomial order Up, i.e. Wp (x Hp. 

In the destriping method, we interpolate 'signal-only' data from the template map m, to remove large-scale (low |u|) modes 
before filtering, or equivalently removing modes inside the circle in the center panel of Figure [15] The smallest mode removed 
is L„„„ — the width of the Fourier transform of the smoothing kernel, .^(u), used to generate the template signal map m, . With 
the low-|u| modes removed, instead of digging a trench along the entire vertical band in the Fourier Plane, as in the usual case of 
polynomial filtering, we filter the same width trench Wp, but only for modes with |u| > Lmin- Residuals from filtering bright point 
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Fig. 15. — Demonstration of Fourier plane inteipretation of destriping process using FDS model 8 simulations. In each column, the 150 GHz / map (top; 
colorscale in MJy/sr with saturated pixels within survey region in white) and corresponding 2D auto power spectrum Cu (bottom) are shown for part of the fourth 
quadrant data. Note that in the power spectrum plots, the color scale is the same for each panel, and is logarithmic. Left: DC-level+slope filtered map mo using 
a signal-only simulation. In the power spectrum, vertical bands are an unavoidable consequence of our scan strategy. Center: same as Left, for a signal+noise 
simulation. The l/f noise, which is added in the scan direction, is obvious as the increase in power in the vertical band. Right: Equivalent plots for the destriped 
map and Cu — there is a large reduction in the 1// noise in the destriped map. and no obvious residuals have been introduced due to the choice of and L„„„. 
In the power spectrum plot, we show the interpretation of the trench dug by the polynomial filter of order n,, on the timestream; L,,,,,, indicates the width of K, i.e. 
the angular scale below which modes are preserved. Modes inside the circle are subtracted from the timestream, as in Equation |Al| 

sources (which dominate the signal at large u) can be avoided by masking the brightest sources during the polynomial filtering in 



Equation A2 



The lower right panel of Figure [T5] demonstrates the result when applying the destriping method with a 6th order polynomial 
— filtering has removed much of the noise along the vertical axis. The corresponding map (top right panel) shows that 1 // noise 
has been heavily suppressed, without the introduction of obvious residuals. We note that some of the 1 // noise has been filtered 
inside L,„,„; this is due to the fact that L,„i„ is not in fact a hard boundary in Fourier space, but rather the width of K, which in the 
present case is a gaussian function. Smoothing therefore allows filtering of different modes with u-dependent 'weighting', with 
the weight equal to K evaluated at each u. 

The plots indicate that though the algorithm is implemented entirely in map/timestream space, its effect is equivalent to a 
Fourier plane filter. Tests on simulated data demonstrate that if uniform scan weighting is used, the algorithm is precisely linear 
in nature. 



B. EFFECT OF FILTERING AND MAP PROCESSING ON DIFFUSE PROPERTIES 

In the following sections, the recovery of diffuse properties is tested with signal-only simulations, binned into the same coarse- 
resolution pixels used for analysis of the QUaD data. When comparing input to output values. Equation [4] is minimized as with 
the real data, except that no error bars are present because signal-only simulations are used. All uncertainties quoted in this 
section therefore reflect the scatter introduced by processing effects, and are quantified by the intrinsic scatter term of Equation]?] 

B . 1 . Total Intensity Spectral Index 

To test the recovery of the spectral index, a signal-only simulation of the FDS dust model 8 is generated, with the model 
evaluated at the central frequency of each QUaD band as quoted in Section ]43| The simulated timestreams are passed through 
the mapmaking code in exactly the same way as the real data, thus incorporating the effects of field-differencing and filtering. 

Note that the m2 mapmaking stage requires source detection and masking, but signal-to-noise thresholds cannot be used to 
determine the location of bright sources with signal only simulations. To circumvent this problem, signal plus noise simulations 
of the same simulated input maps are generated, with the corresponding m\ maps used to determine the locations of sources to 
be masked. The resulting source catalog is then used in the m2 mapmaking stage of the signal-only simulations to reproduce the 
effects of the filtering strategy on signal-only data. 

The input average spectral index was found to be a = 3.5, with a recovered value of a = 3.52 and intrinsic scatter CTq = 0.26, with 
the scatter introduced by the data processing and mapmaking steps. The left panel of Figure ]T6] shows the input and recovered 
pixel values of /150 against /loo, while the right panel shows the input and recovered a as a function of b. Both panels show that 
the input spectral index is recovered to within the scatter introduced by the data processing, with a small systematic defecit of 
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Fig. 16. — Signal-only simulations of spectral index recovery. Left: Scatter plot showing input pixel values (black) and those after filtering, field-differencing, 
and destriping (gray crosses). The input and output mean spectral index a is shown above the plot, with the intrinsic scatter introduced by processing effects. 
Right: Spectral index as a function of galactic latitude b. Mean input is shown as black points, recovered is gray, with error bars indicating the intrinsic scatter in 
each h bin. 

< 1 % in a. 



B.2. Polarization 

The effects of field-differencing and filtering on the recovery of the polarization fraction are investigated by using two signal- 
only simulations, again with FDS dust model 8 as the input total intensity, but using assumed polarization fractions of 2 and 5%, 
with the polarization signal purely +Q (in galactic coordinates). From these maps, we simulate QUaD observations and pass the 
resulting timestream through the mapmaking algorithm as in[B| 



<Q1>,,.=2.1+0.0; <Q/I>y_,=5.C 



-0.2+0.0; <U/I>_=0. I±0.0 



<0/I>^ , =2.0±0.0; <Q/I>,„ =5. 1±0.0 



:-0.0±0.0; <U/I>. 
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Fig. 17. — Signal only simulations of galactic-longitude averaged polarization fraction, as a function of galactic latitude. From left to right: 100 GHz Qjl and 
U /I, 150 GHz Q/I and U /I. Crosses (dots) correspond to assumed 5 (2)% polarization fraction simulations. All input polarized signal is in +Q in galactic 
coordinates, with £/ = everywhere. The average recovered values of the polarization fraction are shown above each plot, with the 'error' representing the 
intrinsic scatter introduced by data processing and mapmaking; in the plots of polarization fraction in bins of b, this quantity is shown as an error bar. 



Figure [Tt] demonstrates that for both 2% and 5% simulations, the input polarization fraction is recovered without introducing 
any strong systematic bias. Fluctuations in polarization fraction due to filtering and map processing effects are 0.1-0.2%, and 
thus introduce a small error on the average polarization fraction measured by QUaD, though such systematic effects do influence 
the recovery of polarization angle. 

Figure [18] shows the recovery of polarization angle from the same simulations. Similarly to polarization fraction, the input 
polarization angle is recovered to within Icr of the introduced scatter, both as an average, and within individual bins of galactic 
latitude. From the average over all pixels, the fluctuation in due to processing and filtering the simulated data at 100 GHz is 
17.1° for 5% polarization fraction, or 26.5° for 2% polarization fraction. At 150 GHz, we find 3.2° (13.1°) for 5% (2%). The 
systematic shift of recovered polarization angle is ^ 5° (^ 2°) for the 2% polarization fraction simulations at 100 (150 GHz), 
and ~ 0.2° (~ 0.4°) for the 5% polarization fraction simulations at at 100 (150 GHz). We conservatively adopt the larger of 
these quantities when estimating systematic errors, i.e. the systematic error on the average recovered (f) is 5° at 100 GHz and 2° 
at 150 GHz. 



C. CELESTIAL COORDINATE MAPS 



Figures 19 21 show the celestial coordinate nii maps. 
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Fig. 1 8. — Signal-only simulations of polarization angle <p\ The left two panels are 100 GHz 4> vs b, and the probability distribution of the pixel values as a 
function of 4>. In the former, the error bars indicate the level of intrinsic scatter introduced by processing effects. Gray crosses (black dots) are for the 5 (2)% 
polarization fraction simulations. The numbers above the plots are the average polarization angles in each b bin, with the 'error' indicating the average intrinsic 
scatter due to data processing. In the probability plots, the gray vertical line shows the input (f> = 0°; gray broken (black solid) lines are for the 5 (2)% polarization 
fraction simulations. The right two panels are the same for 150 GHz. 
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